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ABSTRACT 

We observed the embedded, young 8-10 M Q star AFGL 490 at subarcsecond resolution with the Plateau de 
Bure Interferometer in the C n O (2-1) transition and found convincing evidence that AFGL 490 is surrounded 
by a rotating disk. Using two-dimensional modeling of the physical and chemical disk structure coupled to 
line radiative transfer, we constrain its basic parameters. We obtain a relatively high disk mass of 1 Mg and a 
radius of ~ 1 500 AU. A plausible explanation for the apparent asymmetry of the disk morphology is given. 
Subject headings: circumstellar matter — line: profiles — radiative transfer — planetary systems: protoplane- 
tary disks — stars: formation — stars: individual: AFGL 490 



1. INTRODUCTION 

The study of very young massive stars is of partic- 
ular importance for star formation since it is not yet 
firmly established by what proce ss, disk accretion or stel- 
lar merging, such stars form (e.g ., lYorke & Sonnhalter 2002; 
iDobbs. Bonnell. & Clarkl 120051) Previous authors have 
demonstrated that massive stars can be formed via accretion, 
leading to the formation of large (< 10000 AU) and massive 
(< 1 Mq) circumstellar disks. However, most of these stars 
are located more than 2 kpc away, which makes it challeng- 
ing to confirm unambiguously the presence of accretion disks 
around single objects with radio interferometry due to the lack 
of spatial resolution. 

We focus on AFGL 490 - one of the best-studied young 
stars that are in a transition stage to Herbig Be stars. This 
very young star is located nearby (1 ±0.3 kpc. ISnell et alJ 
QUI, has a mass of 8-10 M©, a bolometric luminos- 
ity of - 2 x 10 3 L Q (B2-3 star, [Panagial fl973T) . is still 
embedded in its parent mo lecular cloud (Ay ~ 40 mag, 
lAlonso-Costa & Kwanll9 89). and drives a high-velocity out- 
flow {Mitchell et al. 1995). Interferometric measurements by 
e.g. Schrever et aLj_j|2 QQ2|. hereafter Paper I) as well a s IR 
observations of Campbell. Persson. & McGregor ( 1986) and 
iBunn. Hoare. & Drew! ( 119 95) have revealed that the star is lo- 
cated inside an ionized ^100 AU cavity cleared up by a fast 
stellar wind and surrounded by a < 500 AU disk-like struc- 
ture, which is enshrouded in a 22 000 AU x 6 000 AU en- 
velope. All these structures are further embedded in an even 
more extended envelope. 

In this Letter, we present subarcsecond-resolution PdBI ob- 
servations of AFGL 490 that prove the existence of a rotating 
disk around AFGL 490 and determine its orientation, size, and 
mass by using comprehensive physical, chemical, and line ra- 
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diative transfer modeling. 

2. OBSERVATIONS WITH THE PLATEAU DE BURE 
INTERFEROMETER 

We observed AFGL 490 in C n O (2-1) at 224.714389 GHz 
using the PdBI in its B configuration in December 2003 (base- 
lines 61.5m-330.5m). The phase reference center of our mea- 
surements was a 20 oo = 03 h 27 m 38 s 55, 5 20 oo = +58° 46' 59".8, 
and the source velocity was set to V\ SI = -13.4 km s" 1 . For the 
line measurements, we used one correlator unit with a total 
bandwidth of 20 MHz and 512 channels (velocity resolution 
of 0.052 km s" 1 ). Two other correlator units with a bandwidth 
of 320 MHz each were utilized to measure the continuum. 

The bandpass and phase calibration were performed on 
the objects 3C454.3 and 2145+067. The additional calibra- 
tions of the phase and amplitude were obtained by observ- 
ing the objects 0355+508 and 0224+671 every 20 minutes. 
For the final phase calibration and data reduction the Greno- 
ble Software environment GILDAS was applied. The maps of 
5 1 2 x 5 1 2 square pixels with . 1 " pixel size were produced by 
the Fourier transform of the calibrated visibilities using nat- 
ural weighting. The synthesized beam size is 0.89" x 0.77" 
(= 890 AU x 770 AU at 1 kpc) with a positional angle of 73°. 
For reduction of the continuum, we excluded the line emission 
from the data. Finally, the continuum was subtracted from the 
original data in the Mv-plane, as it was detected only at the 
source position. 

3. RESULTS 

3.1. Continuum measurements 

Within the 22" primary PdBI beam at 1.3mm, AFGL 490 
was detected point source with a total flux of 1.41 Jy (peak 
value is 0.68 Jy beam" 1 ). The continuum is dominated by the 
thermal dust emission because at 1.3mm the free-free radia- 
tio n contributes less than 5% to the total flux (ps 69 ± 3 mJy, 
see lCampbell etalll9 86). 
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3.2. C O (2-1) line measurements 

The integrated C n O (2-1) line intensity map overlaid with 
the 1.3mm continuum is shown in Fig.^. The continuum 
is barely resolved at the 5a level within the < 2000 AU area 
around the star, and has a beam-like circular shape. In con- 
trast, the line intensity map reveals a ~ 4000 AU arc of emis- 
sion that is centered on the continuum peak and has a PA of ~ 
15°. Furthermore, the C n O emission above the 50% intensity 
level appears as a m 700 AU x 2000 AU bar-like structure 
with the main axis shifted by 400 AU to the northeast from 
the center. 

The red- and blue-shifted parts of the integrated C 17 
2-1) line intensity overlaid with an //-band speckle image 
Alvarez et al. 2004) and the 1.3mm continuum is presented 
in Fig. Both lobes have their maximum 700 AU away from 
the 1.3mm continuum point source, and the line connecting 
them is shifted 260 AU to the northeast (PA of the axis per- 
pendicular to the connecting line is w 26°). Due to the simul- 
taneous measurement of both line and continuum in the same 
bandpass, the positional offset is certainly real. Remarkably, 
the position of the red- and blue-shifted C 17 (2-1) peaks 
perfectly coincides with the position of the CS (2-1) peaks re- 
ported in Paper I. The PA of the C n O emission is smaller than 
the ~ 45° positional angle of the 20000 AU bar-like structure 
(see Paper I), but the emi ssion structure is or thogonal to the 
large-scale CO outflow (Mit chell et alJfT995l) and coincides 
with the PA obtained from the infrared pol arized map a nd the 
2cm radio emission (Campbell et al. 1986; Haas et al. 1992). 
Note that the C n O emission pattern (Fig.[lJ)) is characteristic 
of outflows or rotating disk-like configurations. To distinguish 
between these two possibilities, we analyze the corresponding 
position-velocity diagram below. 

3.2.1. Disk mass 

We estimate the mass of the AFGL 490 "inner disk" 
from the 1.3mm co ntinuum emission, using Eq. 2 from 
iHenning et ail (120001) . We assume the standard gas-to-dust 
mass ratio of 100, a mean disk temperature of 100 K, and a 
dust mass absorption coefficient n c j n at 1.3mm between 1.99 
and 5.86 cm 2 g _1 (grains without icy mantl es, gas de nsities 
10 6 a nd 10 8 cm" 3 , respectively, see lOssenkonf & Henning l 
The resulting H2 gas mass is between 0.8 and 
2.3 M Q and, thus, the source-averaged H2 column density is 
/V(H 2 ) (1 . 3lnm) =l-20-10 23 cm- 2 . 

The gas mass can be also derived from the C n O data. Due 
to missing flux of ^30 % in the vicinity of AFGL 490, even 
in the reverse map side of the red- and the blue-shifted line 
emission, the total integrated flux value is smaller than the 
sum of individual contributions from the red- and blue-shifted 
emission. The sidelobes of the synthesized beam are not con- 
tributing to the low integrated fluxes due to their low intensi- 
ties (<20 %). Thus, we focus on the sum of the red- and the 
blue-shifted emission. 

Assuming optically thin C 17 line emission, the gas mass 
can be derived by 

ln T(r ex + 0.9)exp(16.18/r ex ) 
M(H 2 )[M Q ] = 1.5810- 10 - ; V ' 1 exJ (1) 

(l-exp(-r)) 

x X(C 17 0) (D[kpc])2 / S„Av [Jykms" 1 ] , 

J line 

where t = 0.01 is the optical depth of the C 17 (2-1) line, 
T ex = 100 K is the excitation temperature, X(C n O) = 2.5 ■ 10 7 



is the H2 abundance relative to the amount of C n O molecules 
(H2/COIO 4 , 16 0/ 17 0= 2500), and D =1 kpc is the distance 
to the source. The masses of the red- and blue-shifted gas are 
1 .0 and 0.6 M Q , respectively, giving a total value of « 1 .5 M Q 
which matches the estimate obtained above. 

3.2.2. Disk velocity field 

The gas velocity distribution across the red- and blue- 
shifted C n O emission (solid line in Fig.^) is presented in the 
position-velocity (PV) diagram in Fig.|2] The PV map shows 
a number of single gas clumps, indicating an inhomogeneous 
disk structure. This distribution is not caused by optical depth 
effects as demonstrated by radiative transfer calculations (see 
Sect. 3.3). To fit the PV diagram, we adopt a simple model 
of a Ke plerian disk with a mass that linearly increases with 
radius ( Vo gel et aiIT9 85) and a stellar mass of 8 ± 1 M . We 
assume a disk radius of 1 600 AU based on the spatial ex- 
tent of the C n O emission and a disk mass between 0.5 and 
2 M . With this model, we fit the inner clumps (R < ±0".7) 
of the the PV map and obtain the best-fit inclination angle 
of 35° ±5° (Fig. |2 solid curve). The fit to the most intense 
clumps in the PV map at R = ±0".7 is not as good. It gives 
i = 25° ±5° and super-Keplerian gas velocities in the inner 
disk. In contrast to massive disks around active galactic nu- 
clei that show a modified Keplerian law with an exponent of ~ 
-0.35 (e.g. JKondratko et alJl2005l) . the corresponding curve 
for AFGL 490 (Fig.|2] dotted curve) would only fit to the outer 
disk part. However, it does not fit the inner clumps. Further- 
more, the four most intense clumps in the PV map around 
-12.6 kms" 1 can be best fitted by a rotating gas ring with a 
radius of 700 AU and similar i = 35° ±5° (dashed line). Such 
a small inclination angle is in contrast to the apparent C 17 
morphology implying an inclination angle larger than 60° for 
any circularly-symmetric configuration, e.g. a disk (see el- 
lipse in Fig.Q^). However rotation curves for any disk or ring 
model with larger inclination angles, i > 40°, can only be well 
fitted assuming a low-mass central star, M + < 2 M Q which is 
definitely not the case for AFGL 490. Thus, we propose that 
there is a massive, clumpy disk around AFGL 490 with nearly 
face-on orientation (i ~ 30°) and Keplerian rotation. 

3.3. Disk modeling 

The PdBI C n O (2-1) spectral map of the AFGL 490 disk 
is presented in Fig. [3] We fit thes e spectra u sing the "step- 
by-step" modeling approach ( Semenov et al. 2005, for details 
see Semenov et al. 2006, in preparation). Briefly, we app ly 
the 2D flared-disk model of Dulle mond & DominiM (12004) to 
simulate the disk physical structure, a gas-grain chemical net- 
work to calculate time-dependent abundanc es, and a 2D l ine 
radiative transfer code of Pavlvuchenkov & Shustovl (120041) to 
synthesize the C n O spectra. We assume that the disk has 
the power-law surface density E(r) = Eo('"o / r) p and the ve- 
locity profile V(r) = V()(ro/r) s . The dust grains have MRN- 
like size distribution and the gas-to-dust mass ratio is 100. 
Moreover, we assu me an age of ~0.1 Myr due to the lack of 
reliable estimates. IChurchweiH {l999) reported the dynami- 
cal outflow time of 1.8 10 4 yr which is an lower limit (see 
iHenning et alJ20 00i). The best fit to the observed C 17 (2-1) 
lines is obtained with the following disk parameters: (1) incli- 
nation and positional angles of w 30°, (2) inner and outer radii 
of < 400 AU and 1 400 AU, respectively, (3) a surface density 
gradient p ~ -1.5, and (4) a Keplerian-like velocity profile, 
s pa -0.5. The best-fit disk mass ranges between ~ 0.2 and 
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1 M Q . Note that the modeled line profiles are indeed opti- 
cally thin, t < 0.01. All these values are in a good agreement 
with the observational results, and the values derived from the 
PV diagram. 

4. DISCUSSION AND CONCLUSIONS 

Why do we see such a clumpy disk structure? This ob- 
ject seems to be relatively young and could still be in a per- 
turbed state that remains from an earlier, non-steady accretion 
phase. The viscous dissipation timescale for the AFGL 490 
disk is estimated to be ~ 1 Myr (Pringle 1981), which far ex- 
ceeds its age. This hypothesis is further supported by the fact 
that the ~ 10 4 years old, large-scale CO outflow consists of 
single moving gas clumps, and thus accretion is indeed non- 
steady (Mitchel l et aljfl995h . Another explanation could be 
that the C 17 (2-1) emission traces the densest parts of the 
spiral arms that are easily excited in massive disks by grav- 
itational instabilities. Indeed, the morphology of the C n O 
emission resembles the density structure of a self-gravitating 
disk w ith doubl e arms (upper corner in Fig.^) calculated by 
iFrqmang et alJll2 004al). Moreover, the nume rical simulations 
by Fromang, Balbus, & De Villiers (2004b) imply that such 
disks should also develop a dual structure composed of an in- 



ner thin Keplerian disk fed by a thicker self-gravitating torus 
with nearly uniform rotation. Finally, a clumpy disk struc- 
ture could be a result of a recent encounter with a nearby star 
or due to the gravitational interaction with a wide low-mass 
companion. The complex multiple outflow systems seen in 
the close vicinity of AFGL 490 strongly support this idea (see 
Fig. 4b in Paper I). 

In this Letter, we present clear evidence that ~ 10 M Q stars 
can be surrounded by Kepler ian disks in their earliest evo - 
lutionally phase (see also She pherd. Claussen & Kurtzl200ll) . 
In contrast to the previously reported huge disks around such 
stars with radii of < 10000 AU (e.g., Chini et al. 2004, but 
see Sako et al. 2005) the AFGL 490 disk is much smaller, 
R ~ 1 500 AU. Using advanced theoretical modeling, we con- 
strain basic disk parameters: (1) the inclination and positional 
angles are ~ 30°, (2) the surface density profile is ~ — 1.5, (3) 
the mass is ~ 1.5 M Q (« 50% uncertainty), and (4) the disk 
rotation is close to Kepler's law. 



We acknowledge the help of the IRAM staff both of the 
Plateau de Bure and Grenoble. 
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FIG. 1 . — (a) Map of the total integrated C 17 (2-1) intensity (thin lines) overlaid with the 1.3mm continuum (gray image with dotted lines). The contours 
of the integrated C 17 emission cover intensity levels from 20% (1.5<r rms) to 90% of the peak value (0.25 Jykms~' beam -1 ) in 10% steps; the 50% contour is 
shown by a thick line. The contour lines of the continuum emission correspond to 20% (3a rms), 50%, and 90% of the peak value (0.6 Jy beam ). (Top left) 
This chart shows the midplane density distribution of a self-gravitating massive disk iFromang et al. 2004a, Fig. 9) inclined to 35°. The ellipse indicates the 
sky area covered by such an inclined disk, (b) Integrated red- and blue-shifted C 17 line emission overlaid with the lowest contour as shown in the left panel. 
For the blue-shifted lobe, solid thick contours correspond to gas with velocities between —15.5 and —13.4 kms -1 , while for the red-shifted lobe, dashed-thick 
contours indicate gas with velocities between —12.5 and —9.5 km s~' . Levels are 30, 50, 70, & 90% of the emission peaks (blue: 0.09 Jy beam -1 km s~' ; red: 
0.12 Jy beam -1 kms 1 V The gray area is the //-band speckle image from Alvarez et al. 1 2004). The straight solid line indicates the cut for the position- velocity 
diagram shown in Fig.121 The dashed lines show the major and the minor axes of the emission distributions in the sky plane, whereby the minor axis with different 
position angles (PA) is assumed to be in line with the outflow axis. 
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FIG. 2. — The position-velocity diagram of the C l7 (2-1) emission along the straight line shown in Fig.Qx The thick solid line represents the Keplerian disk 
model with Mdisk = 1 Mq and ( = 35°, while the dashed line corresponds to a model of a 700 AU gas ring with the same orientation. The thick dotted curve 
indicates the disk model with non-Keplerian rotation, V (r) oc r~ 035 . 



